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FLIGHT ASSESSMENT OF AN ATMOSPHERIC TURBULENCE MEASUREMENT SYSTEM 


WITH EMPHASIS ON LONG WAVELENGTHS 

Richard H. Rhyne 
Langley Research Center 


SUMMARY 

A flight assessment has been made of a system for measuring the three com- 
ponents of atmospheric turbulence in the frequency range associated with air- 
plane motions (0 to approximately 0.5 Hz). Results of the assessment indicate 
acceptable accuracy of the resulting time histories and power spectra . Small 
residual errors at the airplane short-period and Dutch-roll frequencies (0.5 
and 0.25, Hz, respectively), as determined from in-flight maneuvers in smooth 
air, would not be detectable on the power spectra. However, errors of ±0.3 to 
±0.6 m/sec (±1 to ±2 ft/sec) at approximately 0.25 Hz can be present in the 
time history of the lateral turbulence component, particularly at the higher 
altitudes where airplane yawing motions are large. An assessment of the quan- 
tities comprising the vertical turbulence component leads to the conclusion 
that the vertical component is essentially accurate to zero frequency. 

Calculations based upon the performance after each flight of the inertial 
platform used to measure low-frequency airplane motions indicate that maximum 
possible linear trend errors on the order of 1.5 to 1.8 m/sec (5 to 6 ft/sec) 
in 10 min could be present in the time histories of the horizontal turbulence 
components. A computer experiment indicated that trend errors of this magni- 
tude would affect the power spectral estimate at the lowest frequency (15 240-m 
(50 000-ft) wavelength) to an appreciable extent only if the turbulence inten- 
sity is very light (i.e., turbulence standard deviation (0 to 10 Hz) is about 
0.9 m/sec (3 ft/sec)). 


INTRODUCTION 

Power spectral methods for designing airplanes have now become widely 
accepted. (See refs. 1 and 2.) As larger and more flexible aircraft fly 
faster, the need for an adequate mathematical model for describing atmospheric 
turbulence at long wavelengths has become more important since the primary 
aircraft-response spatial frequency (cycles/m) varies inversely with airspeed. 

In addition, accurate measurement of atmospheric turbulence over a wide range 
of wavelength, particularly the long wavelengths, is needed to develop a better 
understanding of basic atmospheric phenomena that generate turbulence. (See 
ref. 3.) For these reasons, the Langley Research Center undertook a flight mea- 
surement program referred to as Measurement of Atmospheric Turbulence (MAT). 

The prime objective of the program was to extend power spectra to wavelengths of 
at least 9100 m (30 000 ft) for a variety of meteorological conditions and a 



wide range of altitude, with, the same instrumentation system employed throughout 
the investigation. (See ref. 4 for further details.) 

Measurement of true gust velocities imposes stringent requirements upon 
the instrumentation, particularly when the frequency (and wavelength) range is 
extended over several decades. It is necessary to resolve small fluctuations in 
flow angles and airspeed with a high degree of accuracy. It is also necessary 
to determine accurately the inertial motion of the aircraft as it responds to 
both pilot control inputs and the turbulence field. Since gust velocities are 
determined from the difference between the flow field measurements and aircraft 
inertial measurements (which are related), it is necessary that the instrumenta- 
tion system measure parameters with a negligible phase error. Because of the 
complex nature of the measurement problem, accuracy can best be assessed by 
means of flight-test calibrations. 

The primary purpose of this report is to present the methods used and 
results obtained from a flight assessment of the MAT instrumentation and data 
reduction procedure. This involved taking measurements during pitch-, yaw-, 
and speed-change maneuvers in smooth air (no turbulence) . Special components 
were developed for the instrumentation system. An inertial navigation system 
was modified to enable accurate measurement of the three components of gust 
velocity to a very low frequency. A specially developed airspeed system is 
sensitive enough for measurements of longitudinal gust velocity up to altitudes 
of 19 800 m (65 000 ft). (See ref. 5 for further details.) 

Low-frequency errors (linear trend type), which cannot be assessed directly 
by flight-test calibrations, are discussed, and their effect on power spectra is 
indicated. An artificially generated time history having a known spectral shape 
is used to illustrate the effect on the spectrum of introducing a known trend 
into the time history. The maximum trend errors possible in the horizontal 
velocity measurements from the inertial platform are estimated from postflight 
platform performance. The effect of these errors on spectra of true gust veloc- 
ity is assessed for a range of turbulence intensity. A somewhat similar trend 
error arises from platform measurements of, vertical airplane velocity. A data 
reduction technique is described which removes the trend error from this veloc- 
ity measurement. The technique is based upon independent pressure altitude mea- 
urements and provides essentially trend-error-free time histories of vertical 
gust velocity. 


SYMBOLS 


Measurements and calculations were made in U.S. Customary Units and are pre- 
sented in both the International System of Units (SI) and U.S. Customary Units. 

h pressure-derived altitude based on standard atmosphere table (i.e., 

^ geopotential altitude as defined in ref. 6), m (ft) 

k slope of erroneous linear trend in time history of airplane vertical 

velocity, also zero error in vertical accelerometer which when 
integrated produces this trend, m/sec^ (ft/sec^) 
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scale of turbulence, m (ft) 


horizontal distance between inertial-platform accelerometers and flow- 
direction sensors, m (ft) 

free-stream static pressure, Pa (psi) 

free-stream total pressure. Pa (psi) 

impact pressure, p^ - p, Pa (psi) 

total duration of run, sec 

time , sec 

longitudinal component of gust velocity, positive in direction of 
flight path, m/sec (ft/sec) 

true airspeed, m/sec (ft/sec) 

east-west component of incremental horizontal airplane velocity 
obtained from inertial platform, with arbitrary zero at instant 
that data switch is turned on, positive toward east, m/sec (ft/sec) 

north-south component of incremental horizontal airplane velocity 
obtained from inertial platform, with arbitrary zero at instant 
that data switch is turned on, positive toward north, m/sec (ft/sec) 

incremental vertical airplane velocity obtained from computer integra- 
tion of output from vertically oriented accelerometer mounted on 
inertial-platform stabilized element, positive up, m/sec (ft/sec) 

horizontal airplane velocity along inertial-platform east-west axis, 
positive toward east, m/sec (ft/sec) 

horizontal airplane velocity along inertial-platform north-south 
axis, positive toward north, m/sec (ft/sec) 

lateral component of gust velocity, positive toward right wing, 
m/sec (ft/sec) 

vertical component of gust velocity, positive up, m/sec (ft/sec) 

angle of attack, positive with flow-vane trailing edge up, rad 

angle of sideslip, positive with flow-vane trailing edge toward 
right wing , rad 

incremental vertical acceleration obtained from inertial platform 
with reference to 1g level flight, positive up, m/sec^ (ft/sec^) 



Ah 


P 


AP 

APt 

Alp 


e 


incremental pressure-derived altitude with reference to value at 
beginning of run (see definition of hp) , positive when altitude 
increases, m (ft) 

incremental free-stream static pressure with reference to value at 
beginning of run, Pa (psi) 

incremental free-stream total pressure with reference to value at 
beginning of run. Pa (psi) 

incremental impact pressure with reference to value at beginning of 
run. Pa (psi) 

incremental sensitive airplane heading with arbitrary zero at instant 
that data switch is turned on, measured in horizontal plane, posi- 
tive with nose right, rad 

pitch attitude, measured in vertical plane, positive with nose up, 
rad 


pitch rate measured by body-mounted pitch-rate transducer, positive 
with nose going up, rad/ sec 

wavelength, or distance per cycle, m (ft) 

standard deviation, m/sec (ft/sec) 

roll attitude of airplane with reference to horizontal,’ positive with 
right wing down , rad 

airplane heading, measured in a horizontal plane clockwise from grid 
north, ^ always positive, deg or rad 

yaw rate measured by. body-mounted yaw-rate transducer, positive with 
nose going right, rad/sec 

A bar over a symbol indicates average over the entire run. 


A caret over a symbol indicates that the quantity is given with respect 
to the mean for the entire run; that is, the mean has been subtracted. 


”^Grid north is true north at the platform alinement location, but as the 
platform moves east or west from its initial alinement point , its north-south 
axis is not torqued to point at true north but remains parallel to a vertical ■ 
plane through the meridian at which it was alined. (The north-south and east- 
west axes are torqued to be perpendicular to the local vertical at all times, 
however.) For all practical purposes the inertial— platform axis system can be 
assumed to be alined with true north, considering the latitudes of operation 
and the east-west distances flown in the MAT project. 
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BACKGROUND 


The method of measuring true gust velocities from a flow-vane instrumented 
aircraft was developed simultaneously by both NACA and Cornell Aeronautical 
Laboratory in 1954 (refs. 7 and 8). The velocities so obtained were the true 
motions of the air particles relative to an Earth-fixed axis system. The axis 
system was generally alined with the mean flight path of the aircraft , with gust 
velocity components perpendicular to the mean path of the vehicle obtained in a 
vertical plane and, in later investigations, in a horizontal plane. Power spec- 
tral procedures for the calculation of airplane response to turbulence had been 
developed prior to this time (ref. 9), and the newly developed method of measur- 
ing true turbulence input found ready application. Results of the measurements 
of true gust velocities were immediately shown to be most useful in the form of 
power spectra, as indicated in references 7 and 8. 


Although power spectral measurements of true gust velocity from a variety 
of meteorological conditions and for range of altitude were' desirable , extensive 
recording and combining of the rather large number of corrective airplane motion 
measurements required for data processing was impractical at the time. The 
method was well established, however, and its basics are still in use today. 


As digital methods for recording large numbers of data channels on magnetic 
tape were developed and modern high-speed computers for data processing became 
available, additional power spectral turbulence measurement programs were under- 
taken. (See refs. 10 to 12, for example.) Results from these programs have, in 
general , led to the conclusion that the Von Karman mathematical model of atmo- 
spheric turbulence (described in ref. 13) is valid for the slope of the power 
spectrum (-5/3) at spatial frequencies (inverse wavelength) above about 
3.3 X 10“3 cycle/m (10"3 cycle/ft) for a number of meteorological conditions. 
Considerable uncertainty still exists as to the appropriate value of L, the so- 
called integral scale value, and even to the applicability of the Von Karman 
model at lower values of spatial frequency (longer wavelengths). Experimental 
evidence at longer wavelengths has been lacking because drift-free airplane 
attitude and velocity measurements and long sample lengths required to obtain 
statistical reliability of the power spectral estimates (which must be made at 
smaller frequency intervals) were difficult to obtain. Also selection of appro- 
priate filtering and processing techniques is required . 


CALCULATION OF GUST VELOCITY COMPONENTS 

Time histories of the three gust velocity components were determined in the 
data reduction process from the following equations, which also indicate the mea- 
surements required : 

Longitudinal component 


Ug = (-V) + (V^j^ sin ip + Vgy cos ip) (1) 
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Lateral component 


/s 


(V3) + (-V Aijj + cos ^ - V^y sin \1) + jt^lJ + VoKt) 


Vertical component 


(2) 


A\ 

Wg = (Va) + (-V9 + - ve<)>) 


(3) 


The first term on the right side of each equation represents airflow measure- 
ments. Subsequent terms represent airplane motion corrections. The equations 
are basically those given and discussed in reference 13 » with modifications to 
incorporate the inertial platform measurements directly. (These modifications, 
along with a more detailed derivation of the equations can be found in ref. 14.) 
The assumptions of reference 13 nnd the present study are that (1) all distur- 
bances are small enough to allow use of the angle in radians in place of the 
sine or tangent of the angle, (2) structural flexure between the inertial- 
platform measuring station and the airflow measurement location on the nose boom 
is negligible,' and (3) induced upwash effects on the flow— vane indications are 
negligible . 

The horizontal axes for equations (1) to (3) B.re fixed by the mean heading 
angle "ij/ for each separate run. In calculation of each gust velocity compo- 
nent, fluctuations are given with respect to an arbitrary mean value. In the 
case of the longitudinal component u , for example, V is true airspeed with 
respect to the mean true airspeed for the run, is incremental inertial 

velocity along ^ the platform's east-west axis with mean removed, and V^y is 
incremental inertial velocity along the platform s north-south axis with mean 
removed . The only airplane motion measurements which do not come from the iner- 
tial platform are IJ) and 9, which are obtained from body— mounted yaw— rate and 
pitch-rate transducers, respectively. The distance H associated with these 
measurements is the distance between the platform-mounted accelerometers and the 
flow-direction sensors. The terms involving I account for the difference in 
linear velocity between the platform measuring station and the a and 9 mea- 
suring station due to yaw rate and pitch rate, and are generally quite small. 

The misalinement between the airplane body axis and the horizontal platform axis 
is thus ignored in this case. However, the Va<() term in equation (2) and the 
Vg(t) term in equation (3) are cross terms which do account for the angular mis- 
alinement between the two axis systems with regard to attitude. These second- 
order terms are quite small and can be ignored in most cases . 


INSTRUMENTATION AND HARDWARE 

Since the primary purpose of the MAT project is to obtain accurate measure- 
ments of gust velocity to very low frequencies (long wavelengths), an inertial 
navigation system modified for this application measures airplane motions. (See 
ref. 5 for details of instrumentation.) The airplane motions are then applied 
as corrections to the airflow measurements to obtain gust velocity. 
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The airflow measurements are made on a nose boom extending well ahead of 
the airplane to minimize the effects of the airplane flow field . The nose boom 
can be seen in figure 1, a photograph of the MAT project airplane. Figure 2 is 
a close-up photograph of the special pitot-static head, with balsa flow-direction 
vanes installed. The flow vanes are located 2.01 m (6.59 ft) ahead of the origi- 
nal nose, or 8.37 m (27.45 ft) ahead of the wing leading edge. To insure that 
deflections in the structure between the flow vanes and the inertial platform 
cannot cause angular and/or vibratory effects which cannot be accounted for , it 
is desirable to have a nose boom whose natural frequency is 10 to 12 Hz or 
greater . The tapered construction of the present installation resulted in a very 
stiff, lightweight nose boom, with a natural frequency of about 25 Hz when mea- 
sured on the airplane with pitot-static head and flow vanes installed . 

Details of the individual measurement channels, including the associated 
electronic equipment, method of recording, frequency response, range, and esti- 
mated accuracies are given in reference 5- The sensitive airspeed system devel- 
oped for use in the MAT project is described herein because of its importance in 
obtaining the longitudinal turbulence component. The principle of operation is 
first to obtain the steady-state true airspeed at the beginning of each data run 
by means of conventional measurements of impact and static pressure , and then to 
measure incremental pressures by means of separate auxiliary transducers during 
the remainder of the run. These auxiliary transducers then must cover only the 
range of fluctuations caused by the turbulence and by the pilot in controlling 
the airplane, and thus the required pressure range is considerably reduced, with 
a resulting increase in accuracy and resolution. The method of operation is 
illustrated by the schematic drawing in figure 3* Valves A and B are ordinarily 
open during climb and speed changing so that changing pressures are equalized on 
both sides of the sensitive transducers (shown below valves A and B) and in the 
chambers shown. At the beginning of a data run, the pilot's data switch acti- 
vates solenoids which close valves A and B to "lock" reference pressures which 
exist at that instant in the volume chambers. These volume chambers are heavily 
insulated and thermostatically controlled so that thermal drifts do not cause 
detectable pressure changes on the sensitive transducers during the data- 
recording interval . Once valves A and B are closed , incremental total pressure 
and incremental static pressure are recorded for the remainder of the run in 
addition to coarse-resolution impact pressure and free-stream static pressure . 

Although the pilot maintains a reasonably horizontal flight path during the 
run, small changes in altitude would contaminate the measurements of Ap^ (and 
thus the measurements of incremental true airspeed) at a low frequency if not 
removed in the data reduction by use of the incremental static pressure record- 
ing . In other words, what is needed for obtaining true airspeed fluctuations is 
actually Aq^, or incremental impact pressure (equal to Ap^ - Ap) , which is 
readily obtained in the postflight data reduction. The dynamic requirements of 
the static pressure measurements are not stringent since the airplane does not 
change altitude very rapidly. In order to improve the dynamics of the total 
pressure measurement, interchangeable restrictor orifices are installed at the 
front of the boom according to the anticipated altitude range for the particular 
mission. Their effect is to compensate for changes in damping (of the air column 
plumbing system) due to air temperature changes. 
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Although not indicated in figure 3, there are actually two pressure trans- 
ducers for each of the two incremental pressure measurements, one having a 
recording range of ±1720 Pa (±0.25 psi), and the other a range of ±690 Pa 
(±0.1 psi). The more sensitive of the transducers was allowed to go off scale 
at the lower altitudes where its sensitivity was not needed , but it remained on 
scale during use at the higher altitudes. After the first few flights, it was 
found that the pilot even had difficulty in keeping the ±1720-Pa (±0.25-psi) 
range aP*- transducer on scale in the rougher turbulence at low altitude. (As 
the airplane climbs because of a slight nose— up condition, the airspeed tends 
to decrease; and the increasing altitude and the decreasing speed couple to 
reduce the total pressure; the opposite effect, of course, takes place for a 
slight nose-down out-of-trim condition.) This problem was solved by combining 
the APt and Ap transducer signals electrically so that they were subtracted 
onboard, rather than in the postflight data reduction, to obtain . The 

±1720-Pa (±0.25-psi) recording range for Aq^ then proved to be adequate since 
the effect of the altitude change was no longer present in this quantity. It 
should be pointed out that the linear range of the individual Ap{. and Ap 
transducers was more than ±3450 Pa (±0.5 psi), so that a transducer linear range 
exceedance was not expected to occur . As a precaution , ap^ was recorded and 
monitored over a range of ±3450 Pa (±0.5 psi). In addition, the ap pressures 
were recorded separately, as before, for use in obtaining fine-resolution alti- 
tude changes during the runs. 


ASSESSMENT OF MAT SYSTEM BY IN-FLIGHT MANEUVERS 

Pilot-induced maneuvers in smooth air provide a good practical means of 
checking the overall instrumentation and the data reduction procedure. The 
quantities recorded during the maneuvers are simply inserted into equations (1), 
(2), and (3). Time histories of the resulting turbulence components Ug, Vg, 
and w are expected to remain nearly zero during the maneuver. The maneuvers, 
of course, must be made in very calm air or in a sufficiently large uniform mov- 
ing air mass so that appreciable changes in air movement do not occur during the 
maneuver. In practice, exactly zero gust velocities are seldom attained because 
of the stringent nature of the test . It should be pointed out that the relative 
amplitude and phase of the various motions induced by the pilot are not exactly 
the same as those caused by actual turbulence . For example , to generate appre- 
ciable vertical airplane accelerations and associated vertical airplane veloci- 
ties (the V term) in smooth air maneuvers, the pilot must induce by means of 
elevator control much larger pitch attitude changes (the Ve term) than are 
present during turbulence flight. In addition, the Ve oscillation caused by 
the maneuver strongly correlates with the V„^ oscillation; this is not the case 
for flight in turbulence . The technique employed by the pilot in rough air can 
also influence the magnitude of some motions to a considerable extent . The 
pilot ^s instructions for rough air flight are (1) to maintain airplane attitudes 
reasonably centered within the available recording range by means of gradual con- 
trol movements (special pilot displays connected to the instrumentation provided 
assistance in this task) and (2) to maintain altitude and airspeed as close as 
practical to the altitude and airspeed at the beginning of the run. No special 
attempt is made by the pilot to control the short-period and Dutch-roll motions 
of the airplane . 
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Pitching Oscillation 


Figure 4 presents the result of a typical pitching oscillation maneuver at 
an altitude of 7990 m (26 200 ft). The terms in equation (3) for w are plot- 
ted separately with the same units and scale factor so that their rexative ampli- 
tude and phase can be easily observed. The frequency of the motions is about 
0.5 Hz, which approximates the airplane's short-period stability mode which is 
excited in rough air. The amplitude of the 0.5-Hz pitch attitude oscillation 
(the Vt term) is considerably greater here than that experienced in rough air 
at 0.5 Hz, although the total range covered (about 17.4 m/sec (57 ft/sec)), 
including the gradual nose-down change in trim, is not too unusual for a 10-min 
turbulence run. The total range of in rough air is mainly a function of 

the pilot's skill and attention in controlling the airplane. The iq term is 
also several times larger than that usually experienced in rough air. The 
computed vertical gust velocity w is shown at the bottom of figure 4 and 

should of course be zero. Instead of the ±0.6 to ±0.9 m/sec (±2 to ±3 ft/sec) 

shown . 

The residual gust velocity shown in figure 4 is not considered excessive 
in view of the rather large amplitude of the induced pitching motions. It is 
seen that the lower frequency component of the maneuver (as evidenced by the 
downward trend of the Ve oscillations) is completely counteracted by the 
downward trend of as can be seen by the absence of any low-frequency trend 

in Wg. 

At least one-half of the residual gust velocity oscillation can be accounted 
for by the upwash created at the flow-vane measuring station by the flow field of 
the airplane. Physically, the flow field around the oscillating wing extends up- 
stream far enough to cause the vane measuring incremental angle-of-attack o to 
read high. Calculations for the MAT airplane under average flight conditions, 

based on the method of reference 15, give an induced upwash factor of about 

10 percent of the angle of attack. Close inspection of figure 4 indicates that 
the residual gust velocity is nearly in phase with Va (w appears to be lag- 
ging slightly) , and that a 10-percent amplitude reduction in a would reduce 
the gust velocity oscillation considerably. However, this simple correction 
will not suffice for turbulence measurements, since it is based on quasi-static 
flow effects about the wing. Penetration of the turbulence flow field by the 
airplane, the dynamics of flow buildup about the wing due to turbulence, and 
propagation forward to the flow vane make such an upwash correction considerably 
more complicated. Such corrections have generally been ignored in the past. 

Some justification for not attempting such a correction here is the fact that 
a "noise hump" at 0.5 Hz has not been discernible on the power spectra obtained 
to date with the MAT instrumentation . The lack of such a power peak can probably 
be attributed to the relatively low amplitude of the short-period motions of the 
MAT airplane. (Generally, the amplitude of pitching motions at 0.5 Hz is not 
greater than about one-third that shown in fig. 4.) 


Yawing Oscillation 

Figure 5 presents results from a yaw oscillation made in smooth air at an 
altitude of about 4660 m (15 300 ft). The terms of equation (2) are plotted in 
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the same way as in figure 4. The pilot forced the first cycle of the oscillation 
by rudder movement and then allowed the motion to decay naturally. The rather 
slow decay rate is evidence of the airplane s lack of damping in ynw, which is 
even more pronounced at higher altitudes. The amplitude of the various terms is 
quite similar to that caused by actual turbulence. It might be noted that the 
relative size of the two platform horizontal velocity terms (i.e., cos 'll 

and V sin iJ) during an overall heading change (which occurs here) depends upon 
the airplane ^s mean heading 'll. In the example of figure 5, it^is apparent that 
the aii^lane is on a predominantly north— south heading , since changes 

with ^ with practically no change occurring in Vgy sin 'll. On a 45 mean head- 
ing, both terms would be affected equally. 

The results of the yawing oscillation are more important than the pitch 
results, since airplane yawing motions of approximately this amplitude are con- 
stantly present in turbulence of moderate intensity. In figure 5, lateral gust 
velocity v is ±0.3 to ±0.6 m/sec (±1 to ±2 ft/sec). This erroneous oscilla- 
tion at a frequency of about 0.25 Hz (approximately the Dutch-roll frequency of 
the airplane) will probably be discernible in the time history of lateral gust 
velcity. The size of a resulting hump in the power spectrum at 0.25 Hz can be 
estimated by assuming that the noise is a sine wave with peak amplitude of 
±0.6 m/sec (±2 ft/sec). The mean square, or power spectral area contribution, 
would thus be 


'sine wave 


( Amplitude 


When the amplitude is ±0.6 m/sec (±2 ft/sec), the mean square would be 
0.18 m2/sec^ (2 Vt^/sec^). If moderate intensity turbulence is assumed to have 
a standard deviation a of 3 m/sec (10 ft/sec) , the total area under the power 
spectrum would be 9 m^/sec^ (100 ft^/sec^). The percentage area contribu- 

tion of the noise arising from the yawing oscillation would thus be only 2 per- 
cent and would probably be obscured by normal fluctuations of the power 
estimates . 

The reason for the ±0.3 to ±0.6 m/sec (±1 to ±2 ft/sec) residual gust 
velocity obtained in the yawing maneuver can probably be attributed to a small 
phase difference between the A'ji snd g time histories which has not been 
accounted for. It is noted that the overall heading change (caused by the 
pilot's first rudder oscillation being slightly unsymmetrical ) is compensated 
for quite adequately by the V cos ip term, so that zero gust velocity is 
effectively maintained at this lower frequency. 


Speed Change 

On one or two occasions a speed-change maneuver was performed as a matter 
of interest and to serve as a direct check on the true airspeed measuring capa- 
bility. (Pitching and yawing maneuvers were generally performed in smooth air 
on every flight and were considered an adequate overall check of the instrumen- 
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tation.) In figure 6 the terms of equation (1) are plotted during a speed-change 
maneuver^ Since it was made on an almost exact east-west heading, only the 
^ax ’I' term is active. A very small amount of turbulence is indicated during 
this run by the meager high-frequency content of the incremental true airspeed 
(V term), and of the resulting gust velocity component u The amplitude of 

the speed-change maneuver could not have been much larger at this altitude with- 
out exceeding the ±1 720-Pa (±0.25-psi) recording range of Aq^. 

The results of the speed-change maneuver shown in figure 6 primarily serve 
to show that the special measurements and data reduction procedures for obtaining 
incremental true airspeed are correct. Such a maneuver is not itself representa- 
tive of motions encountered in turbulence . Possible error in the longitudinal 
gust component which correlates with the maneuver is not apparent ; however , very 
small errors could be masked by the very low- intensity turbulence present . The 
steep inertial speed change (V sin ijj) was caused by the pilot applying an 
abrupt power change . Changes of this steepness do not ordinarily occur during 
the turbulence runs since power settings are changed only slightly, if at all. 
Changes of this nature in true airspeed V might occur because of wind shear or 
mountain wave effects, however. Observation of the motion records in turbulence 
indicates that the airplane itself does not respond appreciably to air motions in 
the longitudinal direction until somewhat longer wavelengths are reached . 


ASSESSMENT OF VERY LOW-FREQUENCY AND TREND ERRORS 

Although higher frequency system errors can be assessed by means of in-flight 
maneuvers, extremely low-frequency drifts, or trend errors, cannot be easily 
determined from such tests. The tests would need to be of at least 10-min dura- 
tion and would thus cover a distance of 110 km (60 n. mi.) or more. Even on very 
calm days, the atmosphere is in continual motion, and one probably could not 
expect to travel that far without encountering changing air motions of a minimum 
of 1.5 to 3-0 m/sec (5 to 10 ft/sec). It is apparent, therefore, that low- 
frequency trend errors cannot be readily distinguished unless they are of a sys- 
tematic nature or of appreciable magnitude. 

Almost all of the individual measurements from the MAT system were ade- 
quately assessed for very low-frequency and trend errors by ground-based pro- 
cedures. For example, the flow-direction vanes were locked in a fixed position 
by their calibration jigs, and their electronic outputs observed for stability 
over long periods of time . Likewise the pressure measurements were observed for 
possible trends caused by thermal drifts and other instabilities. Three mea- 
surements, the airplane velocities obtained from the inertial platform, can only 
be assessed for low-frequency errors in an indirect manner. These measurements 
are very important in determining the low-frequency part of the turbulence power 
spectra, and are discussed in the following sections. 


Vertical Airplane Velocity 

Since inertial platform systems are not generally designed to provide navi- 
gation information in a vertical plane, special vertical velocity signals 
extracted from such systems may not be reliable to zero frequency, and conse- 
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quently, they can have undesirable trend errors. The basic difficulty is that 
the vertically oriented platform-mounted accelerometer (from which the vertical 
velocity information must be derived) is susceptible to changing Coriolis forces 
from run to run, as a result of changing ground tracks with respect ^to the 
Earth's rotation. This effect is not compensated for by the system's computer. 
Another effect is the change in the Earth's gravity constant with altitude. 

This effect is largely removed by a special compensation provided in the MAT 
instrumentation which generates an electrical signal proportional to altitude 
from a static pressure measurement . The result of either or both of these 
effects is an apparent shift in the zero of the vertical accelerometer from run 
to run, which when integrated produces an unknown linear trend in the derived 
vertical airplane velocity, and, in turn, in the vertical turbulence component. 
(The accelerometer zero for any one run, although not precisely known, is quite 
stable since the pilot maintains the heading and altitude within fairly close 
limits . ) 

One way of overcoming this difficulty is to integrate the accelerometer 
output in the postflight data reduction rather than to utilize onboard electronic 
integration. Points along the time history where the vertical airplane velocity 
is estimated to be zero are chosen to start and stop the integration. These 
start and stop times are chosen by inspection of the time history of the sensi- 
tive incremental static pressure Ap (an altitude indication); that is, points 
are chosen where the slope of the time history is zero . A mean for the acceler- 
ometer output for this time range is then determined on the computer. When the 
mean (which is determined to a large number of digits and not limited to the 
resolution of the individual acceleration readings) is subtracted and the inte- 
gration performed , the resulting velocity time history has thus been forced to 
start and stop at zero . 

This procedure has limitations however. It is not easily automated into 
the data reduction process, and exact points along the time history where zero 
vertical velocity occurs are hard to determine because of a certain amount of 
noise present in the time history of the sensitive static pressure. (This rela- 
tively high-frequency noise results from cross flow over the static port loca- 
tions on the pitot-static head, particularly in severe turbulence.) Another dis- 
advantage is that zero vertical velocity points may not occur at convenient loca- 
tions near each end of the time history, and, in fact, for the shorter pitch and 
yaw maneuvers, they may not occur at all. 

A more efficient procedure , which accomplishes the same result but with 
greater accuracy and no limitation on the starting and stopping points, is now 
described . The basis for the procedure is the fact that when an erroneous trend 
of slope k is integrated over total time T (with mean removed), a parabolic 
error with a maximum value at T/2 of -kT^/S results. The maximum error at 
T/2 is determined experimentally in the present case by double integration of 
the platform-mounted accelerometer output to obtain a time history of inertial 
displacement. This displacement is then compared with the time history of the 
pressure-derived altitude. An illustrative sketch of the procedure is as 
follows : 
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It should be noted that the end point of the time history of ff Aa^Cdt)^ has 
been forced to agree with the time history of Ah at the end of the run by an 
adjustment to the initial condition of the first integration; that is, the last 
integration is then /(V^^ “ ^az C)dt, where C = Ah (T)/T. The whole pro- 
cedure has been automated in the data reduction process. The displacement error 
at the midpoint of the run is obtained by averaging over ±1/2 sec to minimize 
possible effects of the previously mentioned high-frequency noise on the time 
history of sensitive pressure altitude. The value of k obtained is applied 
as a detrend slope correction to the vertical airplane velocity obtained from 
the first integration. 


As a check on the overall procedure, the vertical airplane velocity 
was again integrated after the detrending procedure and was compared with the 
pressure altitude. The resulting time histories, which are for a 12-min turbu- 
lence run at an altitude of 13 100 ra (il3 000 ft), are shown in figure 7. This 
particular turbulence run is of interest because of the unusually large altitude 
excursions. The amplitude of the excursions obtained by the two methods agrees 
quite well . Certainly all systematic parabolic error between the two quantities 
has been eliminated. Thus, the airplane vertical velocity measurements must be 
essentially accurate to zero frequency. The cause of the very slight discrep- 
ancy between the two quantities near the beginning of the run is not known. 

Other measurements indicate that a large wind shear (decrease in head wind) was 
encountered at about that time and caused the pilot ^s indicated airspeed to drop. 
Such a drop in airspeed might have caused him to put the airplane's nose down, 
causing the 300-m (1000-ft) drop in altitude. Altitude oscillations generally of 
smaller amplitude than the first oscillation of figure 7 are quite often present. 
These are the result of the airplane 's phugoid or long-period mode . Some pres- 
sure noise (discussed previously) caused by the turbulence can be seen about 
5 min after the beginning of the run , and again at 8 min . 
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Horizontal Airplane Velocity and 

The accuracy of the fine-resolution platform incremental velocity measure- 
ments V and V is dependent upon the general performance of the inertial 
platform^nd its coJaputer system. As indicated previously, the accuracy of these 
measurements at very low frequencies cannot be readily assessed to the necessary 
precision by in-flight measurements. A good indication of the platform s per- 
formance and of the maximum probable errors in these quantities is obtained by 
monitoring the velocity readouts after each flight while the platform system is 
still operating. The velocities theoretically should be zero at this time, since 
the airplane is motionless. Errors have built up during the flight however, and 
cause a very small ocillation in the stabilized element with a period of 84.4 min. 
This 84.4-min oscillation is the result of what is referred to as Schuler tuning 
of the platform-computer system. (See ref. 16.) The amplitude of this Schuler 
oscillation is not damped and generally does not decrease during a flight but 
continues to build up with time as errors occur in the system (ref. 16). The 
platform oscillation is too small to be detectable in the pitch and roll mea- ^ 
surements, but is observable in the velocity measurements because of the Earth s 
gravity vector affecting the platform-mounted north-south and east-west 
accelerometers . 

Since the amplitude of the inertial platform's Schuler oscillation is judged 
to be a good indication of the maximum possible error which could have occurred 
during that particular flight, the Schuler oscillation was monitored after each 
flight by making 84-min recordings of velocity. These velocity time histories 
were almost perfect sine waves (of 84-min period) but were sometimes offset from 
zero by an overall trend or "ramp effect." This ramp effect would only be of 
concern with regard to the coarse-resolution absolute velocities and ^ Vy 

and never was of sufficient slope to produce a noticeable effect on the fine- 
resolution incremental velocities used to obtain the turbulence time histories. 
(The ramp effect was generally never as large as the Schuler oscillation, even 
when the inertial system had been in its flight mode for 4 or 5 hr.) 

The magnitude of the postflight Schuler- oscillation for all the flights was 
generally between 3 snd 5 m/ sec (10 and 15 ft/sec) from peak to peak. For the 
purpose of estimating the maximum possible errors in ^ ^^d V^y during tur- 

bulence runs (and therefore in the turbulence time histories), it is assumed that 
the in-flight errors are sine waves having the 84.4-min period_and the postflight 
Schuler amplitude. It is also assumed that the mean heading ^ is such that the 
Schuler error is reflected in the gust component in the worst possible manner ; 
that is, that either sin \|) = 1 , or cos i|) = 1 . It is recognized that the in- 
flight errors are not necessarily smooth sine waves and that error in the iner- 
tial system may not build up in a smooth fashion over the duration of the flight. 
It is believed however, that calculations based on these assumptions provide an 
indication of the maximum possible error in any one run. A 10-min turbulence 
run would have the maximum error, or the maximum effect of trend slope , if it 
occurred at the position on the Schuler error cycle shown in the following 
sketch, that is, 5 min on either side of the zero crossing: 
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The trend can be approximated by a straight line, with the amount of "tilt" A 
given by 


A 



sin 



(4) 


The standard deviation of a straight-line trend with zero mean , which 

has been digitized with closely spaced values, is 


A 

°trend ~ (5) 

V12 

(This is used subsequently to illustrate the effect of the trend on the 

power spectrum.) The following table summarizes maximum possible errors esti- 
mated by the procedure outlined for various Schuler amplitudes: 


Schuler 

amplitude , 

A in 

10 min 

Standard 

deviation 

r (zero 

to peak) 

(eq. 

(4)) 

of trend , 





'^trend 

(eq. (5)) 

ra/ sec 

ft /sec 

m/sec 

ft/sec 

m/sec 

ft/sec 

1.52 

5.0 

1.11 

3.64 

0.320 

1.05 

^2.29 

^7.5 

®1.66 

^5.46 

®.479 

^1.57 

3.05 

10.0 

2.22 

7.27 

.640 

2.10 

3.81 

12.5 

2.77 

9.09 

.799 

2.62 


^Upper limit of Schuler amplitude for most flights. 


These values are considered the maximum possible errors which could be present 
in the longitudinal and lateral turbulence time histories. 

It has been observed that atmospheric turbulence is sometimes embedded in 
a gradually changing horizontal wind field, so that the time history has a 
trendlike appearance . In such cases it is impossible to determine whether a 
trend of the magnitude shown in the table is "noise" or real data. A simple 
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detrending procedure applied to all time histories of longitudinal and lateral 
gust velocity is therefore judged to be undesirable . 


Trend Effects on Power Spectra 

Since the primary purpose of the MAT project is to determine power spectra 
vrtiich are reliable to very long wavelengths, it was desirable to determine the 
extent to which trend errors would affect the spectra. An artificially generated 
time history available from a previous study (ref. 17) was used for this purpose. 
Linear trends of known magnitude were added to the artificial time history, and 
power spectra obtained . These spectra were then compared with the spectrum of 
the original time history. 

The results of the study are presented in figure 8. The spectrum without 
trend is a simulated Dryden transverse spectrum (see ref. 13) of atmospheric 
turbulence, with a value of L of about 180 m (600 ft). The horizontal scales 
have been converted to cycles/m (cycles/ft) (reciprocal of wavelength) by assum- 
ing a 190-m/sec (375-knot) true airspeed to approximate the average true airspeed 
of the sampling airplane. The artificial data were processed by a Fast Fourier 
Transform (FFT) version of the Blackman-Tukey algorithm, the method presently 
being used to process the flight data. The spectral window width (at the 1/2 
power point) is 0.024 Hz and is approximately that required for processing 10-min 
data runs and maintaining 30 statistical degrees of freedom."^ 

The upper curves in figures 8(a), (b), and (c) are the spectral results 
obtained from the turbulence plus trend time histories when the ratios, 

''trpnd/'^ turbulence* I-O* 1-5. respectively. (The appropriate slope 

of the linear trend which was added to the artificial turbulence time history in 
each case was computed by eq. (5).) The spectral results for the turbulence 
alone , that is , for the time history without trend , are shown as the lower curve 
in each figure for comparison . The lowest frequency power estimate shown in 
these logarithmic plots at an inverse wavelength of about 1.64 x 10“^ cycle/m 
(5 X 10"° cycle/ft) or at a wavelength of 60 960 m (200 000 ft) is actually the 
value computed at zero frequency and is generally discarded for actual data . 
Although values at zero frequency have customarily been discarded, they probably 
should be retained in cases where effects of trend error are known to be quite 
small. (The symbols in fig. 8 indicate the power estimates at the five lowest 
frequencies . ) 

For the general upper limit of the possible trend error given in the previ- 
ous table (^tpend = 0.479 m/sec (1.57 ft/sec)), the turbulence standard deviation 


^The magnitude of the random fluctuations observed in the power spectrum 
obtained from the artificial time history is not directly comparable with flight 
data with regard to calculated statistical degrees of freedom. This is because 
the synthetic time history was generated by the addition of a large number of 
sine waves, whose amplitudes were adjusted to provide the desired- spectral shape 
The magnitude of the random fluctuations could be reduced to an arbitrarily low 
level by simply increasing the number of sine waves or by increasing the length 
of the time history. 
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would have to be as low as 0.957 m/sec (3-1^ ft/sec) turbulence = 

before the effect could be as severe as that shown in figure 8 (a; . Similarly, 
if the turbulence had a standard deviation of only 0.479 m/sec (1.57 ft/sec) 

(°^tr?nd'^'^ turbulence “ trend effect could be as severe as that depicted 

in figure 8(b;. Figure 8(c) is generally unrealistic but might apply to one or 
two flights for which an extremely large Schuler error was recorded in combina- 
tion with a turbulence run of very low intensity. Generally, turbulence with 
a standard deviation of less than about 1 m/sec (3.5 ft/sec) is considered 
extremely light, and unless it is very consistent, of very long duration, or of 
particular interest for some other reason, it would probably not be considered 
worth evaluating. Therefore, on the basis of figure 8, the trend would affect 
the power spectral estimate at the lowest nonzero frequency (15 240-m (50 000-ft) 
wavelength) to an appreciable extent only if the turbulence intensity is very 
light . 

It is believed that if a narrower spectral window width were employed in the 
analysis, the effect of the trend would be considerably less- on points adjacent 
to zero frequency. (The effect might be greater on the zero-frequency point, but 
as previously mentioned this value would probably be discarded in any event.) In 
the case of actual data, employing a narrower window of course requires that the 
time histories be longer to maintain the same statistical confidence . The use of 
a more rectangular shaped spectral window, available from the so-called frequency 
averaging algorithm (see ref. 18), should also reduce the effect of the trend 
considerably at the 15 240-m (50 000-ft) wavelength. 


CONCLUDING REMARKS 

A flight assessment has been made of a system for measuring the three com- 
ponents of atmospheric turbulence in the frequency range associated with airplane 
motions (0 to approximately 0.5 Hz). Results of the assessment indicate accept- 
able accuracy of resulting time histories and power spectra , 

In-flight maneuvers in smooth air were used to determine residual errors 
in the time histories of gust velocity at the airplane short-period and Dutch- 
roll frequencies of approximately 0.5 and 0.25 Hz, respectively. The small 
errors obtained are not expected to be detectable in the power spectra . A 
residual error at about 0.25 Hz with amplitude of ±0.3 to ±0.6 m/sec (±1 to 
±2 ft/sec) is, however, expected to be detectable in the time histories of lat- 
eral gust velocity, particularly at the higher altitudes where airplane yawing 
motions are large . 

A data reduction procedure which effectively determines a more accurate 
zero for the inertial-platform vertical accelerometer and in turn removes the 
accompanying trend error in the airplane vertical velocity term is described. 

An assessment of the quantities comprising the vertical gust velocity measure- 
ments, including a comparison of the integrated vertical airplane velocity which 
has been trend-corrected, with pressure-derived altitude for a 12-min run, leads 
to the conclusion that the vertical turbulence component is essentially accurate 
to zero frequency. Calculations based upon the inertial-platform performance 
after each flight are used to indicate that maximum possible linear trend errors 
on the order of 1.5 to 1.8 m/ sec (5 to 6 ft/sec) in 10 min could be present in 
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the time histories of the horizontal g-ust velocity components. A computer exper- 
iment indicated that trends of this magnitude would affect the power spectral 
estimate at the lowest frequency (15 240-m (50 000-ft) wavelength) to an appre- 
ciable extent only if the turbulence intensity is very light (i.e., turbulence 
standard deviation (0 to 10 Hz) is about 0.9 m/sec (3 ft/sec)). 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
October 15, 1976 
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Figure T-- Time histories of integral of vertical airplane velocity and pressure-derived 

altitude in a 12-min wind-shear turbulence run. 
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Figure 8.- Linear trend effect on pover spectrum of artificially generated 
time history. (Symbols denote power spectral estimates at five lowest 
frequencies . ) 
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Figure 8.- Continued 
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Figure 8.- Concluded. 
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